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A THEORETICAL AND EXPERIMENTAL INVESTIGATION OF 
CYLINDRICAL ELECTROSTATIC PROBES AT ARBITRARY INCIDENCE 

IN FLOWING PLASMA 

By Renaldo V. Jenkins and W. Linwood Jones, Jr. 

Langley Research Center 

SUMMARY 

This report pi'esents the theory for calculating the current collected by a negatively 
biased cylindrical electrostatic probe at an arbitrary angle of attack in a weakly ionized 
flowing plasma. The theory was constructed by considering both random and directed 
motion simultaneous with dynamic coupling of the flow properties and of the electric 
field of the probe. This direct approach yielded a theory that is more general than static 
plasma theories modified to account for flow. 

Theoretical calculations are compared with experimental electrostatic probe data 
obtained in the free stream of an arc-heated hypersonic wind tunnel. The theoretical 
calculations are based on flow conditions and plasma electron densities measured by an 
independent microwave interferometer technique. In addition, the theory is compared 
with laboratory and satellite data previously published by other investigators. In each, 
case the comparison gives good agreement. 

INTRODUCTION 

Simplicity makes the electrostatic or " Langmuir 1 ' probe attractive for plasma 
diagnostics; however, the theoretical interpretation of its current- voltage characteristics 
is extremely complicated. Nevertheless this problem has been extensively analyzed and 
for the case of collisionless plasmas at rest has been adequately described. Unfortu- 
nately, the theory for interpreting the response of an electrostatic probe in a plasma with 
directed flow is considerably more complicated than the static case. 

For years it was argued that sufficiently long cylindrical probes alined parallel to 
a flowing plasma responded independent of the flow. On a general scale, the credibility 
of this argument was severely weakened by Bettinger and Chen’s publication of L. H. 
Brace’s Langmuir probe data from the Explorer 17 satellite. (See ref. 1.) During this 
experiment, a probe with a length- diameter ratio of 410 exhibited an unexpected increase 
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in current as the axis of the probe was brought into alinement with the satellite's velocity 
vector. (See fig. 1.) Since this result could not be explained by any static plasma 


Current, pA 



Figure 1.- An example of ion current peaking as published by Bettinger and 
Chen in reference 1 (L. H. Brace's observation). (Fit for BP 595 from 
ref. 1 with 10$ H + and 90$ 0 + . The solid line represents the theoretical 
curve; the x's, the experimental curve with 4 - K/r; the dots, the 
experimental curve with (j> = K/r^. ) 

theory, Bettinger and Chen (ref. 1) developed a new theory for flowing plasma. In their 
analysis the ion current peaking was explained as an end effect due to the finite length of 
the probe. Subsequent plasma measurements were performed by Hester and Sonin 
(refs. 2 and 3) under laboratory simulation of the satellite-plasma conditions. In these 
experiments large ion current increases, when the probe was alined with the flow (see 
fig. 2) were also attributed to end effects. 
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Figure 2.- End effect at a probe potential of the order of the stream energy. 

(Fig. 9 of ref. 3.) U^KTg/n^ )'°' 5 = 38 , \/ r ^ = e ^p/ KT e " “ 26 °* 

There were other flowing plasma measurements, however, which exhibited no cur- 
rent peaking and in which a rigorous static theory (ref. 4) modified for use in supersonic 
flow (refs.. 5 to 7) gave convincing agreement with experimental results. 

At the time when the present authors became cognizant of these references (refs. 1 
to 3), they had completed both published (ref. 8) and unpublished work where probe end 
effects were not considered. In each case, the electron density was inferred by using 
both electrostatic probes and a microwave interferometer, and the agreement between 
the two techniques was good. 

The present research was therefore undertaken to resolve whether an end effect 
was present for the conditions of the previous work and also to develop an improved 
generalized theory for calculating the positive ion current collected by a cylindrical 
electrostatic probe oriented at an abritrary angle of attack in a flowing plasma. 

SYMBOLS 

A,B,C empirical coefficients defined in text 

d noncollection length (see appendix A) 

E ratio of bulk flow velocity to ion velocity 

E,E r electric field intensity; radial component 

e electronic charge 


\ 
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empirical factor and coefficient 
empirical factor 
current 

current to an infinite probe 

magnitudes of electron and ion current densities at surface of a probe 
Boltzmann’s constant 
length of probe 
mass 

average electron number density along microwave transmission path 
ion or electron number density (assumed equal) 
magnitude of electronic charge 
probe radius 
temperature 
time 

bulk flow velocity 
voltage 

ion velocity (approximately the speed with which ions enter the side of the 
Debye shield), 

permittivity of free space 

angle of attack between probe and flow direction (called angle of orientation) 
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D 


4>r 
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Debye length, 


fr 0 Kr e y /a 

^ e 2 n / 


angle of noncollection (see appendix B) 
difference of potential between probe and plasma 

dimensionless potential difference between probe and plasma, e(V - Vp^KT ( 
angle of collection (see appendix B) 


Subscripts: 


e,- electron 

end total current contribution at zero orientation angle 

i,+ positive ion 

m maximum 

min minimum 

p plasma 

r retarded 

1,2 transition points 

3 point at which maximum angle of collection has its minimum value 

THEORY 
Other Theories 

A detailed review of theory and application of Langmuir probes can be found in the 
references by Loeb and de Leeuw (refs. 9 and 10). The original analytical expression 
for the ion current collected by an electrostatic probe oriented at an angle 6 to plasma 
flow is the one by Mott-Smith and Langmuir (ref. 11). Reference 3 gives that expression 
as 


5 



1/2 


= 2enUr p L|sin 2 0 - 
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(1) 


This expression, however, is only valid if the ion thermal energy is much less than the 
flow-directed kinetic energy (kt. « (l/2)m.U 2 ). In the case of parallel alinement 
(0 = 0), the expression reduces to 


J co = 2enr p L| 


f-2e<p r 


v l/2 


nn 


( 2 ) 


It so happens that this equation is also the expression for the ion current collected 

by a probe in a static plasma in the orbital motion limit (e<£> p » Id\ » -i m.U 2 j. For 

these conditions an infinite cylindrical probe alined parallel to flowing plasma collects 
ions as if the plasma were stationary. 

Although this expression is inappropriate for the general case of ion collection by 
a probe parallel to a flowing plasma, it led to the application of other static theories to 
flowing plasmas. One of the most significant and realistic of these static theories is the 
one by Laframboise (ref. 4). This theory was developed for a two-component (ions and 
electrons) plasma, where each component had a Maxwellian velocity distribution. 

The theory is complex and need not be thoroughly covered here, but some salient 
points are worth noting. The case of negative probe potential is considered primarily, 
since it furnishes the most information pertinent to the present investigation. 

By using the mathematical framework of Bernstein and Rabinowitz (ref. 12), 
Laframboise (ref. 4) derived expressions for the absolute magnitude of the ion and elec- 
tron current densities at the surface of a probe. For negative potentials these expres- 
sions are for the ion current density 



(3) 


and for the electron current density 


i 


e 


en 



(4) 


From this last expression Laframboise (ref. 4) showed that the electron tempera- 
ture could.be determined from a semilog plot of j e against V; this result agreed with 
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Langmuir's original conclusion. He was able to obtain this result since for slightly 
negative potential 



( 5 ) 


The electron temperature is therefore easily obtained from this theory but the ion 
number density is not. 

From equation (3) 


e 


fZW ill. 

2’ im j IVxp-Tj 




The calculation of n from this equation is complicated by the fact that the dimen- 
sionless current I’ is dependent on n through Thus, n can only be computed 

through a tedious iterative procedure. Fortunately, Sonin (ref. 5) developed an algorithm 
for calculating n which was much more satisfactory than the iterative procedure. 

Sonin's algorithm was evaluated experimentally in supersonic flow (refs. 5 to 7) 
and gave convincing results. The publication of Brace’s observation by Bettinger and 
Chen (ref. 1), however, invalidated its applicability in certain flowing plasmas. Subse- 
quently, Hester and Sonin (ref. 5) also confirmed Brace's results in a wind tunnel. 

In references 1 and 3, theoretical analyses are included which explain the observed 
peaking of the ion current. Although they each attributed the peaking to end effect, they 
disagree on certain points. Furthermore, neither theory is generally applicable to 
flowing plasma. Therefore, the present theory was developed to be applicable in a more 
general sense to flowing plasma. 


Present Theory 

A theory applicable to flowing plasma must account for flow dynamics. In the 
present theory flow dynamics have been included by the consideration of such things as 
flow velocity, probe angle of attack, flow-directed kinetic energy, side collection, angle 
of collection, and so forth. The inclusion of such dynamic considerations has resulted in 
a theory widely applicable to flowing plasma. 

This report restricts the theory presented to the negative potential difference 
(0 < 0) region of the probe response curve. (See fig. 3.) Although this theory should 
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Secondary ionization 



Figure 3*- Complete Langmuir probe characteristic. Ion current exaggerated. 
Theory of present report restricted to region between points A and B. 


be extendable to the region of positive <|>p, such an extension is beyond the scope of this 
work. 


The theory is relegated to two distinct classes. The names and criteria for these 
classes are as follows: 


Class I - Thermal: With probe at an angle of attack of 0°, the probe current must 

be dominated by particles that enter the Debye shield as a result of random thermal 

1 2 

motion. Thus, the particle flow-directed kinetic energy - nijU* is much less than the 
thermal energy KTp 


Class II - Flow directed: With probe at an angle of attack of 0°, the probe current 

must be dominated by particles that enter the Debye shield as a result of flow-directed 

motion. Here the kinetic energy -i rmU^ is much greater than the thermal energy KTj. 

Z 


Notice that the criteria do not exclude particle contribution by thermal or flow- 
directed motion from either class. This nonexclusion is mandated by the physical reality 
that no flowing plasma is without both types of motion. 
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The theory was developed by use of the following steps: 

(1) Assume that the Debye shielding length concept is valid for any charged body 
whose characteristic dimension (probe diameter here) is less than the mean free path of 
the plasma. Thus, only the charged particles within the Debye shield can be affected by 
the probe. 

(2) List the ways in which ions on electrons could enter the Debye shield, that is, 
ions enter the shield's end as a result of flow-directed motion. 

(3) Construct current contribution terms by using the list of step 2, that is, the end 
ion current is enU(r p + X D )*V cos 9 . (The ideal current is the algebraic sum of all such 
terms.) 

(4) Determine whether any of the terms of the ideal current should be eliminated 
for physical reasons. (Example: The flow-directed electrons entering the end of the 
shield have insufficient kinetic energy to overcome the electric field, except for very 
small potential differences. Under this condition there is a physical basis for elimi- 
nating the end electron current contribution terms.) 

(5) Modify the current contribution terms where necessary to account for such 
items as maximum angle of collection, true collection length, and the effect of the various 
velocities on particle collection. 

In some instances the modifications of step (5) took the form of semiempirical 
coefficients and factors. The term semiempirical is used since these coefficients and 
factors are computed from physical properties of the plasma or are the decade normali- 
zation of such computations. 


Thermal (Class I) 

The derived expression for calculating the current collected by a negatively biased 
probe at small angles of attack (0° < 6 < 15°) is 

I = 2en7iYrp + A-d)(L - d + cos 8 ) ( V + ) 

- 2en7r(rp + Xj))(L - d_ r cos 0)(v_)exp(-e<£p/KT e ) 

+ enU(r p + Xj^tt cos 9 - enU^rp^Tr cos 9 

+ 2enL(r p + X D )(sin ^ m ) E (V + )0 F “ 2enL^r p ^U sin 9 (7) 

The first and second terms represent the currents due to ions and electrons, 
respectively, that enter the Debye shield as a result of thermal motion. The effect of 
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flow dynamics is included through the quantities d + and d_ r (ion and electron non- 
collection lengths). These factors must be included because in the time required for the 
charged particles to reach the probe surface, the ions and electrons move downstream 
these lengths. (Expressions for computing and d_ r are derived in appendix A.) 
The third and fourth terms are end currents due, respectively, to flow-directed ions and 
electrons. 


The fifth and sixth terms are the flow-directed currents collected by the probe side 
area. The variables \js m and 6 of the fifth term are introduced to account for the 
modified capture area and efficiency of the probe while at an angle of attack. (See 
appendix B.) The exponent E is used to weight the angle -of -attack effect and is defined 
as the ratio of the bulk flow velocity to the ion velocity, E = -y-. 


It was empirically determined that the coefficient F is related to E by the 
E 

equation F = ^ — where p is the decade in which E occurs. Mathematically, F 

can be regarded as simply the normalization of E, where the normalizing factor is the 
largest E of the decade which contains E. ^Examples: if E = 4.6, then 

F =-£- = M= 0.46; if E = 13.7, then F = -1^1- = 0.685.) 

in 7 ' J 


10 10 


2 ( 10 ) 


The fourth and sixth terms are not present if 


i m e (U sin 9) 2 


|m e (U cos 6) 2 


and 


are each less than e0 p . Usually this is the case; therefore, they can 


be neglected. For large <£ p < 0, the second term can also be neglected. 


At zero angle of attack almost all the probe current is contributed by the first and 
second terms since Class I - Thermal is being considered. Note that the side collection 
(fifth) term dominates the shaping of the ion current as a function of orientation angle 9. 


Flow-Directed Class II Theory 

This class consists of two subclasses determined by the ratio E of the bulk flow 
to the ion velocity. For subclass A this ratio is less than 10, whereas for subclass B it 
is greater than 10. The theory presented for each subclass is applicable in the ion- 
saturation region of the response curve only. (See fig. 3.) 

These subclasses came about as the result of the observation that the empirical 
multiplying factor, in the derived equation for ion current, must be different for the dif- 
ferent E ratio groups. A transition E of 10 was suggested by the fact F (the decade 
normalization E) has a value of 1 for the first time at E = 10. An analytic search for 
further transition E ratios w T as fruitless; thus, it was concluded that only these two sub- 
classes exist. 
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Subclass A.- The ion current is 


I = GenU(r p + A D ) 2 i7 cos 9 + 2FenL(r p + A D )(sin ^ m ) A (sin ’/' min ) B U P (sin 9) C (8) 

Here the effects of flow are included both directly (through U the flow velocity) and 
indirectly (by means of the multiplying factor G). The factor G is an empirically 
determined end collection efficiency defined in a subsequent paragraph. It depends not 
only on flow properties but on probe orientation angle as well. 

As a first cut, one would intuitively surmise the side collection term (term 2) to 
be 2enL(r p + Aq^U cos 9. Note, however, that this form of the term can be true only if 
all the flow-directed ions entering the Debye shield from the side are collected. In 
reality total collection is an impossibility since some of the ions have sufficient momen- 
tum to take them outside the shield. Thus, the angle of collection and the efficiency of 
collection must be accounted for. This accounting is accomplished through the intro- 
duction of the efficiency coefficient F and the sine of the maximum angle of collec- 
tion i// m . As the angle of orientation 6 increases, \f/ m decreases until it reaches 
its minimum value of For orientation angles above the angle where this occurs, 

the side current becomes more directly dependent on 9, and thus sin 6 is included. 

The roles played by the coefficients A, B, and C in the shifting of importance of 
sin ^ m i n , sin 6, and sin i f/ m will become clear later when they are defined. 

The factor G equals F = — for angles of attack between 0° and the angle (in 
degrees) numerically equal to approximately (E/2) + 1. From this angle to the first 
transition point angle 0j, G equals zero. Higher angles of attack (9 > $i) have G 
equal to 1 - F. The first transition point angle is the angle at which the exponen- 
tial coefficients A, B, and C make their first value change. 

The values of A, B, and C are defined as follows: 


A = E ^ 
B = 0 
C = 0 


sin > exp 


log. sin * mln j\ 


E - 1 


The coefficients A, B, and C initially have these values when 0 = 0. 
A = 1 

( log e sin 4/ 



sin \p m < exp( 


min 


V E - 1 


The equations for and ^ min are presented in appendix B. 


min 
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( ( log e sin i// , \ \ 

I sin < exp^ -g— -j j and sin & 2 = sin ^ m J 

Thus it can be seen that there are two possible transition points (from definition (1) 
to (2) and from (2) to (3)). One occurs at 

sin = exp 
and the other at 

sin x {/ = sin 6 2 

(These points are called, respectively, transition points 1 and 2.) Transition point 2 may 
not occur since the conditions that 

sin 0 2 = sin 

and 

sin xj/ < exp 

must be simultaneously satisfied for it to do so. Transition 2 can only take place after 
transition 1 even though the requirement sin xf/ m = sin 0 O is met. 

Ill z 

As is shown in appendix B, each angle of orientation 9 has a maximum angle of 
collection xj/ m associated with it; therefore, the transition point angles can be uniquely 
determined. 

The shaping of the curve of ion current against 9 is determined by the side col- 
lection (second) term. 

Subclass B .- This subclass has an E ratio greater than 10, and the ion current is 
I = [l + (1 - F)]enU(r p + \ d ) 2 ti cos 9 

+ [l + (1 - F)]2enL(r p + X D )(sin ^ m ) A (sin ^ min ) B (sin f?) C (U) F (9) 

The exponential coefficients A, B, and C retain the same definitions as were given 
for subclass A. Here the curve shaping is also dominated by the side collection (second) 
term. 



log e sin xp min 
E - 1 
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EXPERIMENTAL APPARATUS 


Electrostatic Probe Array 

The probe array used in this investigation consisted of four electrostatic probes 
protruding from a 15° half-angle wedge with a rounded vertex angle. Figure 4 gives two 


Top 



(a) Cross-sectional slice of electrostatic probe array. (Side view.) 



(b) Top view of probe array. 

Figure 4.- Electrostatic probe array. 

views of the array. Each probe consisted of two main parts: a 0.0254-cm-diameter 
tungsten electrode, and an insulator holder fashioned from 0.3-cm-diameter glass or 
quartz tubing. In addition, two probes of the array had a third part, an insulator cap. 

The top probe had a 0.127 -cm -diameter quartz cap. (See fig. 4(a).) The adjacent probe 
had a 0.099 -cm -diameter ceramic cap. The capped probes had exposed lengths identical 
to the 1.905-cm exposed lengths of the uncapped probes. Thus, each probe had a length- 
diameter ratio of 75. 
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The two capped probes and one of the bare probes were each fixed biased nega- 
tively with respect to the tunnel walls, and during a run an oscillograph was used to 
record the probe current. This recording was accomplished by using the circuitry of 
figure 5. The remaining uncapped probe was swept biased to record the entire voltage- 
current characteristic. The experimental circuit (fig. 6) utilized a linear voltage sweep 
from approximately -6 V to 2 V during a period of 1 to 10 milliseconds. Both the voltage 


Galvanometer 



Differential 

amplifier 



Figure 6.- Block diagram of swept probe setup. 
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and current waveforms were displayed on a dual-beam scope where they were photo- 
graphed by a high-speed framing camera. 

The number of probes in the array was later reduced to two. Both of these were 
bare, one fixed and the other swept biased. They were located in the two center posi- 
tions of the array. (See fig. 4(a).) The capped probes were removed to eliminate the 
cap shock waves. Visually, only the 0. 127-cm-diameter quartz cap had a detectable 
shock; however, because of the uncertainty in shock interference, both capped probes 
were removed. 


Microwave Interferometer 

A 10-GHz microwave interferometer was used to measure independently the elec- 
tron number density in the plasma flow. A discussion of microwave interferometer 
theory can be found in references 8, 13, and 14; but fundamentally, the interferometer is 
a bridge circuit which compares the phase of the data path signal with that of the refer- 
ence path signal. The data path signal undergoes a negative phase shift relative to the 
undisturbed reference path signal as it transverses the plasma. This phase shift is 
related to the average electron density along the microwave transmission path by the 
following equation: 

N = 7.38 X 10® Ay electrons/cm® (10) 

where Ay is the phase difference in degrees. 

A diagram of the experimental arrangement is shown in figure 7. Note in this 
diagram that the transmitting and receiving horn antennas are surrounded by microwave 
_absor.bing_mater.ial — This material eliminated multipath-propagation-by-absorbing-micro- 
wave energy that could have otherwise been reflected off the tunnel walls. Such multipath 
propagation would have resulted in erroneous phase shift measurement. 

As a further guard against multipath effect, the stainless- steel wedge of the probe 
array was positioned a minimum of 12.065 cm downstream of the microwave beam. This 
wedge position located the probe tips 7.62 cm downstream of the microwave path. 

Rotating and pivoting the probe array at this location did not affect the phase of the 
microwaves in the absence of plasma and thereby indicated that the probe array produced 
no measurable multipath effects. 
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Description of Facility 

A diagram of the Langley 1-foot hypersonic arc tunnel (30.5 cm) is shown in fig- 
ure 8, and a description of the tunnel can be found in references 15 and 16. The electric 
arc -heated air is expanded through a conical nozzle and produces a low- density flow at a 
nominal Mach number of 12 in the test section. An analysis of the accuracy of the test 
conditions in this facility is given in reference 17. 

In this test series, the facility was operated at stagnation pressures of 0.507 to 
0.709 MN/m 2 with stagnation enthalpies of 6.1 to 9.8 MJ/kg. This mode of operation 
resulted in bulk test stream velocities of 2.98 to 3.71 km/sec. 
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Arc heater 



Figure 8.- Schematic drawing of the Langley 1-foot 
hypersonic arc tunnel ( 30.5 cm). 


RESULTS AND DISCUSSION 
Experimental Search for End Effect 

An interesting but inconclusive part of this work was an experimental attempt to 
detect probe end effects. The scheme of detection was based on the hypothesis that in 
the presence of end effect, probes without end current contribution should have curves of 
ion current against angle of attack very different in shape from those with end current. 
Thus, the electrostatic probe array was designed to eliminate the end ion current for 
certain probes. This design was accomplished by capping two probes of the array as 
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previously described. To be effective, these insulator caps had to prevent ions from 
reaching the physical end of the probe and also had to eliminate the ions that would other- 
wise enter the end of the Debye shield. To accomplish this, the diameter of the insulator 
cap should have been a minimum value of twice the sum of the probe radius and the Debye 
length 2(r p + X D ). 

Insulator caps with such diameters were almost certain to produce shock waves 
which could invalidate the experiment. Nevertheless the end blocking criteria were used 
to determine the minimum cap diameter. The 0.099- cm ceramic cap had no visually 
detectable shock; therefore, the experiment proceeded on the assumption that if any 
shock was present, it was a weak conical shock. And as such, it would not appreciably 
change the flow angle. 

During the experiment the angle of attack of the probe array was varied and the ion 
currents for one capless and two capped probes of the array were simultaneously 
recorded. Runs were conducted both with air and with nitrogen as test media. The 
result of one such air test is presented in figure 9. In this typical figure, the ion current 
for each probe is plotted against angle of attack. The diamond, square, and circular 
symbols, respectively, represent data for capless, 0.099-cm-diameter ceramic-capped, 
and 0. 127-cm-diameter quartz-capped probes. The magnitudes of these currents are dif- 
ferent; however, their shapes are strikingly similar. This result seems to be an indica- 
tion of the absence of end effect. This position is further supported by the "valleying" 
present for each of the probe currents. Such valleying is directly opposite to the pro- 
posed end effect peaking of references 1 to 3. 

A comparison of the pertinent test conditions for the present experiments and those 
of references 1 and 3 is given in table I. Theoretical considerations for the cylindrical 


TABLE I.- PERTINENT TEST CONDITIONS 


Source 

n, 

electrons/ 

m 

Neutral-neutral 
mean free path, 
m 

r P’ 

m 

L, 

in 

u, 

m/ sec 

Present 

report 

a 3.9 x 10 15 

[3.5 X 10 15 tv 1,2 x 10 16 ] 

a 5.47xl0‘ 5 

[3.92 x 10" 5 to 5.83 x 10“ ?) 

>5.5 x 10" 4 

1.27 X 10" 4 

1.905 x 10" 2 

a 3.17 x 10 3 

[2.9B x 10 3 to 3.71 x 10 3 ! 

Bettinger 
and Chon 
(ref. 1) 

6.39 x 10 U 

4.107 x 10" 3 

-1236 

2.8 x 10~ 4 

0.23 

7.9 x 10 3 

Hester and 
Sonin 
(ref. 3) 

b (l) 3.19 x 10 11 
b (2) 2.69 x 10 12 
b (3) 2.69 x 10 13 

1.651 x 10" 3 

All mean free 
paths were 
larger than 
the test cham- 
ber dimensions 

1.27 x 10" 4 

5.715 x 10~ 2 

b (l) 7 753 
b (2) 21 232 
b (3) 67 143 


a Value corresixmds to data of figure 9. Bracketed values give range for present test scries. 

^These values not given, but computed by assuming ion temperatures of (1) 2 K, (2) 15 K, and (3) 150 K and T g = Tj x 10^. 
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probe geometry are well satisfied (length to probe radii range from 150 to 821). The 
neutral -neutral mean free paths are several orders of magnitude greater for refer- 
ences 1 and 3 than for the present experiments. The plasma electron densities were 
several orders of magnitude lower in references 1 and 3 as compared with the present 
work. The ratios of probe radii to Debye length were less than unity for references 1 
and 3 and approximately one and a half decades less than the ratio for the present 
experiments. 



Figure 9- - Variation of ion current with angle of attack 
for capped and capless probes. 
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Measured Electron Concentration and Temperature 


The microwave and electrostatic probe studies were made in the free stream of the 
Langley 1-foot hypersonic arc tunnel (30.5 cm). Electron number densities of 3.5 X 10 
to 1.2 x 10 16 electrons/m 3 were measured by using the nonperturbing microwave inter- 
ferometer. The values of electron concentration obtained by applying Laframboise's 
theory (ref. 4) through the use of Sonin’s algorithm to parallel swept electrostatic probe 
data agreed within 10 percent, and suggested the validity of Laframboise's (Sonin's 
algorithm) theory for the flowing plasma. The semilog plotting method was used to 
determine the electron temperature. This method yielded electron temperatures of 
2100 Kto 3870 K. 


Application of Present Theory 

In this section the general theory of this report is applied to data representing each 
class. In each case, the ion current for a fixed bias probe was computed and compared 
with the measured ion current. The test series by the present authors furnishes the data 

for Class I (thermal). Class II, subclass A data are from the report by Bettinger and 

Chen (ref. 1), and the data for Class n, subclass B is that of Hester and Sonin (ref. 3). 

Class I (thermal) .- The data for one air and one nitrogen test are used. The free- 
stream conditions for the air test were as follows: 

(1) n = 3.68 x 10 15 electrons/m 3 (microwave interferometer) 

(2) T e = 3020 K (semilog plotting method) 

(3) V + = 915m/sec 

(4) U = 3.18 km/sec 

(5) The ions were NO + . 

The bulk flow velocity U was not measured, but was obtained from figure 10. 

This figure is a calibration plot of stagnation enthalpy against the free- stream velocity 
taken from the work of Guy (ref. 18) for the Langley 1-foot hypersonic arc tunnel (30.5 cm) 
with air as the test medium. 

The ion current collected by the bare fixed bias (v = -6 volts; = -7.5 volts) 
probe for various angles of attack is shown in figure 1 1 as circular symbols. The solid 
curve of this figure is the ion current computed from theory (eq. (7)). Since the probe 
was negatively biased into the ion saturation region, only the first and fifth terms of the 
equation are significant. 

As can be seen in figure 11, the theoretically predicted current is in good agree- 
ment with the measured current. 
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Listed below are the free-stream conditions for the nitrogen run. 

(1) n = 4.36 x 10 15 electrons/m 3 (microwave interferometer) 

(2) T e = 2761 K (semilog plotting method) 

(3) V + = 1280 m/sec 

(4) U = 3.06 km/sec 

(5) The ions were N + . 

The ion current measured with a fixed bias (V = -6 volts; <pp = -7.39 volts) bare 
probe is shown in figure 12. The solid curve is the calculated current. The value of 
bulk flow velocity U used is slightly high since it was obtained from figure 10 and is 
for air, and nitrogen is the test medium for the results shown in figure 12. (An unpub- 
lished nitrogen curve produced by Guy gives approximately the same velocities for air 
and nitrogen when the stagnation pressure is the same. Lower stagnation pressure as is 
the case for this nitrogen run gives slightly lower velocity.) 



The theoretical current was therefore expected to be slightly high because of the 
flow velocity error. Figure 12 shows the disagreement to never be greater than 10 per- 
cent. The primary point of disagreement is the shape of the two curves. From equa- 
tion (7) it is apparent that the shaping of the theoretical curve is dominated by 6 ^ 

(fifth term). The quantity F equals U divided by 10 V + , and U is thus the shaping 
factor. The lower (true) flow velocity would give better agreement between the shapes of 
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the theoretical and measured currents. (Data for tests using nitrogen with stagnation 
pressure compatible with fig. 10 gave curve shape agreement like that in fig. 11.) 

Class II, subclass A .- The data of Bettinger and Chen (ref. 1) are presented in fig- 
ure 1 (which is fig. 8 of ref. 1). In addition to the data (r p , T e , etc.) of this figure, 

the probe length L and $ p are given directly as 23 cm and 0 p indirectly as 
m i 2\ 

<Pp = 0.25 — U J. Taking the ions as 0+ gives <£ p as -5.17 volts. Applying equa- 
tions (All) and (A12) to the data results in t. and d, values of 9.085 X 10 -7 sec and 

+ + 

7.177 mm. The data can also be used to compute: 

V + = 1069 m/sec 



Ap = 4.107 mm 
F = ^j = 0.739 
| + 1 = 4.69° 

G = 0.739 (0 < 0 < 4.69°) 

By using equation (B3), the minimum angle of collection i^ min is found to be 
3.65°. The orientation angle 0 3 at which i// min occurs is obtained from equation (B2) 
as 37.58°. By solving equations (B7) and (B13), the angles of attack d\ and d 2 at 
which transitions 1 and 2 occur are found to be, respectively, 27.67° and 31.43°. Since 
the first transition point angle has been found, the factor G can now be specified 
for the complete range of 6 as 

r 0.739 (O<0<4.69°) 

G = (0 (4.69° < 6 g 27.67°) 

0.261 (27.67° < 6 £ 90°) 

With these computations completed, equation (8) can be utilized and the theoretical 
ion current for various angles of attack can be calculated. The current so calculated is 
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presented in figure 13. A comparison of the theoretical curve of this figure with the 
experimental data of figure 1 gives good agreement. The slight disagreement that does 
exist is for angles above 5° or below -5® and should be attributed to the fact that the 
experimental data have not been corrected to account for the difference in p for these 
angles of attack. This difference in <p p is a consequence of the slowly varying swept 
voltage applied to the probe; that is, the experimental data were not obtained with a 
fixed 0p. 



Figure 13- - The theoretical ion current computed by present theory for Bettinger 
and Chen's data. (The circular symbols correspond to the dots of fig. 1 and 
the squares correspond to the X's of the same figure.) Note that all the 
data of figure 1 have not been transferred; however, all the extreme points 
have been transferred. 

Class II, subclass B.- Unfortunately, the data of Hester and Sonin (ref. 3) are pre- 
sented (fig. 2) in nondimensional form, and the electron number density is not given. 
Thus, an exact analysis is precluded. However, sufficient data are presented to support 
an analysis of the curve shape to be expected. Reference 3 complements the data of fig- 
ure 2 with the following data: 

Tp = 0.127 mm 



Tj near zero 
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Thus Xj), L, and T e are found to be, respectively, 1.651 mm, 57.15 mm, and 
= T iA 0 ’ 2 - this point it is necessary to assume a value for Tj and to assume 
T e = Tj^/lO -2 . 

Assume that Tj = 2 K (near zero) yields a T 6( of 200 K Since values for the 
electron temperature T e and Debye length Xp are known, the electron number 
density n can be calculated from the definition of Aj). Such a computation results in 
an n of 3.19 x 10 1 * electrons/m 3 . The ratio E of bulk flow velocity U to the ion 
velocity V + is given as 38, and V + is computed* as 204 m/sec and results in a U 
of 7753 m/sec, By proceeding as in the example for subclass A, t+, d + , * min> e lt 
e 2 , and e 3 are calculated as 6.793 X 10" 7 sec, 5.266 mm, 4.09°, 7.07°, 18.58°, and 
19.67°, respectively. 

o 

Theoretical values of ion current can now be computed by equation (9), and non- 
dimensionalized by using the current at zero angle of attack as a normalizing factor. The 
nondimensional ion current thus obtained is presented in figure 14. A comparison of the 
shape of this theoretical curve with the shape of the data curve of figure 2 gives good 
general agreement. An exact agreement was not expected since exact values of Tj 
and T^T e were not known. 



Figure Ik.- The theoretical nondimensional ion current for the data of 
reference 3* = 2 K and T e = (Peak current has a 

dimensional value of 0 . 06>6 microampere.) 


To eliminate the possibility that the agreement for these curves is due to a fortui- 
tous choice of (2 K), two other curves using different values of Tj were computed. 

iThe ions were singularly ionized argon (A + ). 

^The sin should have several decimal places to prevent errors in (sin 
since A can be large; that is, here A = 38, 1, or 0. 
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Assumed ion temperatures Tj of 15 K and 150 K result in electron number den- 
sities of 2.69 x 10 12 and 2.69 x 10 13 electrons/m 3 . With T e = Tj/lO" 2 , the electron 
temperatures are 1500 K and 15 000 K, respectively. The corresponding nondimens ional 
ion current curves are presented in figures 15 and 16. The shapes of these curves are 
also in good general agreement with the data curve of figure 2. 



Figure 3 . 5 .- The theoretical nondime ns ional ion current for the data of 
reference 3. Tj. = 15 K and T e = iJlQ' 2 . (Peak current has a 
dimensional value of I .285 microamperes.) 



Figure l 6 . - The theoretical nondime nsional ion current for the data of 
reference 3, = 150 K and T g - T^/lCT^. (Peak current has a 

dimensional value of 38.534 microamperes.) 
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Electron Temperature Determination Using Present Theory 
Thermal (class I) .- At zero angle of attack 6 = 0°, equation (7) becomes 

I = 2 en 7 r(A D + r p )(L - d + )V + - 2en7r(A D + r p )(L - d_ r )(V_)exp^- 

+ enU(r p + A d ) 2 tt -enU^r p ^ 2 ir 

For the electron-retarding region of the response curve (fig. 3), this equation can 
be written as 


I = 2en7r(A D + r p )(L)V + - 2eiw(A D + r p )(L)(v.)exp(- 
+ 2nU7r^r p + Ap^ 2 - emrU^r p ^ 2 

since neither d + nor d_ r is an appreciable part of L. If 
1+ = 2en*(A D + r p )(L)(V + ) 

I_ = -2eni7(A D + r p )(L)(v_) 

and 

X end = enU *( r p + X D) 2 - enU ( r pf* 

then 

1 = X + + X end + (O exp kt^) 


By taking the natural logarithm and 
A-[log e (l - I + - I end )] = jg-. But since 


Thus 


d log e (l 


and T 


differentiating with respect to 0 p , one gets 
cp = V - V p and V p is constant d<£ p = dV. 

e can be found by plotting the current against 


voltage. On semilog paper, such a plot is a straight line with slope 

KT, 


The electron 


temperature is easily computed once the slope is obtained, e and K being known 
physical constants. 


This plotting procedure is exactly the same as that used by Langmuir, Laframboise, 
Sonin, etc. to find the electron temperature. Even though it was demonstrated only for 
parallel probe alinement, it is a simple matter to show that it holds for arbitrary orien- 
tation angle. Experimentally, the present authors determined the same electron 
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temperatures for various angles of attack. Sonin also obtained the same electron tem- 
peratures for various angles of attack (ref. 6). 


Flow-Directed (Class n) Theory 

The theory for subclasses A and B of Class II is valid only in the ion- saturation 
region of the response curve (fig. 3) where only ions are collected. As a consequence of 
no electron collection, the electron temperature cannot be determined by using the theory 
of Class II. Fortunately, in the electron- retarding region of the response curve, elec- 
trons are again collected. This electron collection produces an electron current hun- 
dreds of times greater than the ion current and means that nearly all the total current is 
due to thermal motion. The criterion for Class I is thus satisfied. This result illus- 
trates how classes can be switched (Class II to Class I only) by varying the voltage, so 
that results in the electron- retarding region of the response curve can be obtained. 
Therefore, the electron temperature Tg can be obtained by the procedure presented 
for class I. 


Class Determination 

The application of the theory for one class to data of another class results in gross 
errors; therefore, the data should be correctly classed with the theory. 

Although final classification must be determined by applying the class criteria at 
an angle of attack of 0°, there are certain other characteristics which are helpful. 

Class I (thermal) current data for flowing plasma plotted against orientation angle d, 
for example, has a valley point at an angle of attack of 0°. A similar plot of Class II 
(subclasses A and B) current data has a peak point at 6 equal zero. The distinctive- 
ness of the valleying and peaking effects is directly related to the bulk flow velocity U 
and the ion velocity V+. Class I data for a slow flowing plasma with high V + would 
have an indistinct valley. In fact, a flow velocity U of zero would result in a straight 
line, since the current is independent of the*>rientation angle for a stationary plasma. 
Higher flow velocities produce more pronounced valleys. The peaking of Class II is 
dominated by the ratio E of the bulk flow velocity to the ion velocity. As the ratio E 
increases, the peaking becomes sharper. 

Recall also that the ratio E is used to determine whether subclass A or B is 
represented. Knowing that certain data are of the Class II type is meaningless by itself. 
The bulk flow velocity U and the ion velocity V + must be known ^or V + is com- 
puted from knowledge of values of T e and m^. With U and V + known, E is 
computed. If E is less than 10, then the data belong to subclass A, and if E is 
greater than 10, they belong to subclass B. 
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APPLICATIONS 


It is readily apparent, from the valleying and peaking characteristics of the curves 
of ion current against angle of attack of this work, that electrostatic probes can be used 
to determine local plasma flow direction. This is done by simply recording the orienta- 
tion of the probe as it is pivoted and rotated until the valley or peak point current is 
reached. With this orientation the probe axis is parallel to the local flow; hence, the 
flow direction is that of the probe axis. 

There is also the possibility of using the theory to determine the ion density and 
flow velocity. 

By using a general curve-fit procedure or a theoretically computed family of 
curves (curves of current against n with E, F, x D , etc. determining the members 
of the family), it should be possible to obtain the ion density. A curve-fit scheme would 
probably be more accurate; undoubtedly, it would be complex and tedious and would 
involve one or more iterative procedures. 

Since the ion current for various angles of attack is highly dependent on the flow 
velocity, one may be able to extract the velocity through some analytical manipulation of 
the ion currents for two or more angles of attack. If n is known, then the coeffi- 
cients E and F can be determined empirically. The definition of E (or F) could 
then be used to obtain the velocity. This second approach, unfortunately, requires an 
independent measure of the ion density, 

CONCLUSIONS 

In the experimental part of this work, a microwave interferometer and an electro- 
static probe array have been used to measure the electron number density in the free 
stream of a 1-foot arc-heated hypersonic wind tunnel (30.5 cm). The electrostatic probe 
technique also yielded electron temperature measurements. Both types of measurements 
were made for a wide range of tunnel free -stream conditions with air and nitrogen as 
test media. 

The analytical part of this work consisted of the development of a general theory 
for calculating the ion current collected by a cylindrical electrostatic probe oriented at 
arbitrary angles of attack in a flowing plasma. This theory is composed of two central 
classes (Class I (thermal) and Class n (flow-directed)). Class II is further divided into 
subclasses A and B. It was not determined whether there are conditions for which both 
Class I and Class II theory are equally valid. Nor were the conditions for switching 
from Class I to Class II obtained. 
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The significant conclusions are as follows: 

1. The present theory is applicable to a wide range of flowing plasma conditions. 

2. By using the present theory, the electron temperature is determined by the 
semilog method previously developed. 

3. The peaking and "valleying" of ion current is not the result of end effect pei se, 
but the consequence of side collection which is highly dependent on angle of attack, and 
the ratio of the flow velocity to the ion velocity. 

4. The theory of this report is in good agreement with Sonin’s algorithm results, 
when applied to the conditions of the test series (parallel alinement only). (No end effect 
was detected.) 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 26, 1974. 
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APPENDIX A 


THE NONCOLLECTION LENGTHS d + AND d 


If a cylindrical probe (see fig. 17), having a potential V, is immersed in plasma of 
effective potential Vp, it will collect charged particles (ions or electrons). The par- 
ticles collected are limited to those within a region extending a Debye length Xj) from 
the surface of the probe at any one time. (The Debye length is the maximum distance at 
which a charged body can affect other charged bodies in a plasma.) Since the surface of 
the probe is cylindrical in shape, the imaginary surface (Debye shield) that encloses this 
region is a cylindrical shell capped by a hemispherical shell. A plasma flowing parallel 
(0 = 0°) to the axis of the probes interprets the Debye shield as a cylindrical shell capped 
by a disk, since it projects the hemisphere as a circular disk. Therefore, this shape is 
assumed for the Debye shield (fig. 17) for all small angles of orientation. 



Noncollection length 


"A 



Figure 17.- Probe with its cylindrical Debye shield. 


The particles outside the Debye shield are not influenced by the probe's presence 
(where it is assumed that the probe creates no shock waves). Upon entering the shield 
they are affected by the probe because its electric field E influences them immediately. 

This electric field E arises as a consequence of the potential difference 
(<pp - V - Vp) between the shield and the probe. Furthermore, it can be concluded that 
the E field only varies radially with respect to the probe. This conclusion is readily 
deduced from the fact that the only possible E field is the one normal to both the equi- 
potential surfaces of the probe and of the shield. The lines of the E field lie along the 
radii joining the probe and the shield, and its strength is 
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APPENDIX A - Continued 


E 

r *D 


(Al) 


The force acting radially on a charged particle that enters the E field is 


F r = qE 


r “ 


q0p 


(A2) 


The acceleration a r that the particles of mass m p experience in the E field 
can be written, by means of Newton's second law, as 


a 


r 


dV _ F r 
dt m p 


(A3) 


which can be rewritten by using equation (A2) as 

dv q^p 

dt ” nipA.j) 

Separating variables and integrating gives velocity as a function of time 


q<M 

v(t) - — f- + c i 

m p X D 


(A4) 


(A 5) 


where is the constant of integration. By applying the boundary condition at time 
zero, the particle enters the shield with initial velocity V 0 and gives the value of the 
constant Cj as V 0 . Equation (A5) can be written as 


V(t) = 


q<V 

+ V 

mpX D 0 


(A6) 


The velocity V(t) can also be written as the time rate of change of the coordinate going 

from the shield to the probe V(t) = — . Making this substitution into equation (A6) gives 

dt 


dr 

dt 


q$ p t 

iyo 


+ v r 


(A7) 


By integrating and applying the boundary condition r = 0 at t = 0 (the instant that 
particles enter the field), it is found that 

,2 


r = 


q^pt 

2m p x D 


+ V 0 t 


(A8) 
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APPENDIX A - Continued 


This equation can be used to calculate the total time t that it takes for the particles to 
travel from the shield to the probe if the distance r is replaced with Xp. If the poten- 
tial difference <p p is negative, then ions are accelerated toward the probe and equa- 
tion (A8) for an ion is 


q<V + 2 

2^^-+V + t + - X D = 0 (A9) 


For an electron, equation (A8) becomes 


-q^-r) 2 

2m e Aj) 


+ V t 


- X D = 0 


(A10) 


By solving equation (A9) for the time (t + ) that it takes ions to travel from the Debye shield 
to the probe surface, the roots are 

t = ~( V +) m i x D x D m i\/ v + 2 + (2q0p/ m i) 

+ q<J> P * q<^p 


The root that is desired is 


m i X D V + X D m j\/ V + 2 + 2q0 p /n^ 

q^p q0^ 


(All) 


The other root is useless since it is negative and negative time has no relevance here. 
Multiplying equation (All) by the bulk flow velocity U gives the length d + for which 
no ions are collected. This length is 


d + = Ut + = 


-VD m i U 

q^n 


X D m 




q<^r 


(A12) 


By solving equation (A10) for the time t_ r it takes a decelerating electron to reach the 
probe, the following roots are obtained: 


*-r- 


m e X D V - ± f m e V - ‘ 

q0 P q<^p V 2 


q<p t 


(A13) 


m e V. 


It is immediately apparent that both roots are imaginary if the kinetic energy 
2 

- of the electron as it enters the shield is less than the potential energy q <£ p of 
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APPENDIX A - Concluded 


the field. Physically, this condition means that no such electron can reach the probe. 

If all electrons entering the shield have kinetic energies less than the potential energy of 
the field, the noncollection length d_ r is exactly equal to the probe length L. (This 
condition requires the second term of eq. (10) to vanish.) Equation (A13) represents two 
real roots when the electron has kinetic energy sufficient to overcome the field and reach 
the probe. However, only the root 


t-r = 


m e X D V - 

<l4>n 


1 


2m X D 

5^— M 


I 

m e V. 


- 


(A14) 


is physically possible. The other root would require the time t_ r to increase as the 
entry velocity of the electron increases. This increase is impossible since less time is 
needed to travel a certain distance as the velocity of travel is increased. 


The noncollection length d_ r for electrons in a plasma of bulk velocity U is 


d„ r = Ut_ r 


Um e X D V - 

<l^p 


N/2m> D U I m p v - 2 ~ 
Q0p 'V 2 “ q *P 


(A15) 


In the special case where <p p = 0, equations (A9) and (A10) reduce, respectively, to 


v + t + -x D = 0 


V_t_ r - X D = 0 


Solving these equations gives t + = Xq/V + an d t_ r = Therefore, when = 0, 

the ion and electron noncollection lengths are, respectively, 


d + = Ut + 


ux D 

V + 


d_ r = Ut_ r 


UX D 

V 
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APPENDIX B 


THE MAXIMUM ANGLE OF COLLECTION \!s 

T m 

A probe oriented so that its axis forms an angle 9 with a plasma flowing with 
velocity U has a component of velocity U sin 6 across it. As a result of the velocity 
component across the probe, not all the ions that enter the Debye shield are collected. 
This result is a consequence of the trajectories that ions in the field must have. Any ion 
entering the shield at a point such as point A of figure 18 would follow a slightly curved 
trajectory. This curved trajectory takes the ion out of the shield where it cannot be 
collected; whereas, if the point of entry is shifted along the shield, from A toward C, a 
point B will be reached where the entry will result in collection. Ions entering between B 
and C all have trajectories which carry them to the probe surface where they are 
collected. 



Figure 18.- Cross section of probe and its Debye shield taken perpendicular 
to the probe axis. AE and. AE' are exaggerated trajectories of 
non collection. 

The angle i^ m formed by the lines joining B and D and C and D is the maximum 
angle of collection. Calculating this angle by using the trajectories of ions is difficult, 
if not impossible. Fortunately, the noncollection angle £ can be calculated easily and 
thus can be found. By referring to the geometry of figure 18, the side of the tri- 

angle opposite the angle £ has a minimum value of Ut + sin 6. The hypotenuse BD of 
this triangle is A D + r p . This angle of noncollection in degrees is 

/Ut sin 6\ 

* = arc (B1) 
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APPENDIX B - Continued 


Since i// m and | must sum to form a right angle, the angle of maximum collection is 


\p m = "° " £ = 90 ° " arc sin 


^Ut + sin 9\ 
. Xd + r p < 


(B2) 


Note that the angle 9 has one unique associated with it and that decreases 

as d increases. This trend is valid until the angle 6 = 9 3 for which \p m equals 
yp . is reached. For all 6 values above 1 9 3 , the maximum angle of collection has 
its minimum value of ^ min . 

By referring to figure 19, the angle can be no smaller than one half the angle 

subtended by the chord created by projecting the diameter of the probe onto the shield. 



Figure 19*- Cross section of probe and shield showing geometry 
used in determining a minimum ralue of 


Therefore, 


i b . = arc sin 
^min 


D + r p ( 


(B3) 


is the minimum value of i// m . Now with ant * equation (B2) the orientation 

angles and 9^, with transitions 1 and 2, can be calculated. Transition 1 takes 
place at 


sin ip = exp - 


/logp sin \p. 


m 


E - 1 


(B4) 
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APPENDIX B - Continued 


By means of equation (B3) 


sin 


exp< 


this equation can be rewritten as 

lo Se[ r p/( x D + r p)]\ 

E - 1 J 


Therefore, 


Af 


m 



(B5) 


(B6) 


Substituting this value of >p m into equation (B2) gives 



90° 


- arc sin 



sin 0, 


+ r v 


(B7) 


Equation (B7) is easily solved for the first transition-point orientation angle Since 

transition 2 occurs when sin = sin 9 2 , then must equal 9 2 and ^m can be 

replaced with d 2 in equation (B2). If 90° is subtracted from both sides of equation (B2) 
after this replacement, then equation (B2) becomes 


o /Ut. sin 0 2 \ 

d 9 - 90° = -arc sin — ± £] 

2 l A D”p J 


Taking the sine of equation (B8) gives 


sin (fl 2 - 90°) = 


-Ut + sin 9 2 
X D + r p 


(B8) 


(B9) 


By using the trigonometric identity, sin(a - b) = sin a cos b - sin b cos a, equation (B9) 
can be reduced to 


-cos 0 2 = 


-Ut + sin 0 2 
X D + r p 


If equation (BIO) is divided by 


sin 0 2 and multiplied by -1, it becomes 


cos e 2 
sin d 2 


= cot (?2 = 


Ut + 

X D + r p 


(BIO) 


(Bll) 


37 



APPENDIX B — Concluded 


Inversion gives 


tan &2 = 


x D +r p 

Ut + 


Thus the second transition- point orientation angle 6 g is 


(B12) 


02 ~ arc tan 


L D 


+ r, 


Ut, 


(B13) 


The orientation angle where the maximum angle of collection reaches its 
minimum value can be found by replacing in equation (B2) by its minimum value 

of equation (B3). 
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